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Radiation from a collapsing object is manifestly unitary 
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The process of gravitational collapse excites the fields propagating in the background geometry 
and gives rise to thermal radiation. We demonstrate by explicit calculations that the density matrix 
corresponding to such radiation actually describes a pure state. While Hawking’s leading order- 
density matrix contains only the diagonal terms, we calculate the off-diagonal correlation terms. 
These correlations start very small, but then grow in time. The cumulative effect is that the 
correlations become comparable to the leading order terms and significantly modify the density 
matrix. While the trace of the Hawking’s density matrix squared decreases during the evolution, 
the trace of the total density matrix squared remains unity at all times and all frequencies. This 
implies that the process of radiation from a collapsing object is unitary. 


Introduction. One of the most pressing problems in 
modern physics is the information loss paradox in black 
hole physics. Since Hawking radiation is purely thermal 
PQ, it is possible to convert a pure state into a mixed 
state, which is forbidden in unitary quantum mechanics 
[2J- It was often argued that subtle correlations between 
the emitted Hawking quanta which are usually neglected 
could be enough to recover information about the initial 
state and convert an apparently maximally mixed ther¬ 
mal state into a pure state This point of view was 

also often criticized by noticing that small corrections to 
the leading order Hawking terms are not enough to re¬ 
cover unitarity [5j. Most of this discussion is given for 
the case of radiation from the pre-existing black hole, or 
in the limit when the horizon is already formed in infi¬ 
nite future. The purpose of this paper is to address this 
problem from a point of view of an asymptotic observer 
who is observing a time-dependent gravitational collapse. 
We find indeed that the process of gravitational collapse 
and subsequent evaporation is manifestly unitary as de¬ 
scribed in the foliation of an asymptotic observer. 

We used the functional Schrodinger formalism which is 
especially convenient for this question since it gives us the 
time evolution of the system rather than radiation from 
a pre-existing black hole ISffiS]. We start with a mas¬ 
sive shell which is collapsing under its own gravitational 
pull. This process induces a non-trivial time-dependent 
metric which then excites the field quanta. The pro¬ 
cess of the gravitational collapse takes infinite time for 
an outside observer, however, radiation is pretty close 
to thermal when the collapsing shell approaches its own 
Schwarzschild radius. Our formalism gives us an explicit 
form of the wavefunction of the emitted radiation, which 
contains complete information not only about the diag¬ 
onal Hawking terms, but also about the non-diagonal 
correlations terms. Correlations between the Hawking 
quanta are at first indeed negligible with respect to the 
diagonal terms. However, time evolution creates progres¬ 
sively more off-diagonal terms than the diagonal ones. 
Moreover, time evolution is such that these cross-terms 
become of the same order of magnitude as the Hawking 
terms. As a result, the density matrix for the emitted 
radiation is significantly modified, in particular it is not 


purely diagonal. We calculate the time evolution of the 
complete density matrix as a function of time and fre¬ 
quency. The relevant quantity that we want to obtain is 
the trace of the density matrix squared (Tr(/5 2 )), which 
tells us whether the system is in a pure or mixed state. 
We find that if we take only diagonal terms in density ma¬ 
trix then Tr(ph 2 ) diverges from unity, which means that 
the state goes from pure to mixed. This is the standard 
Hawking’s result which implies information loss. How¬ 
ever if we include the off-diagonal terms then Tr(p 2 ) re¬ 
mains unity at all frequencies and all times during the 
evolution. This means that the initial state stays pure 
during the evolution. An observer who could measure 
the whole density matrix would not lose any information. 
This is the main result of our analysis. 

The formalism. We consider a thin shell of mat¬ 
ter which collapses under its own gravity. We use 
Schwarzschild coordinates because we are interested in 
the point of view of an observer at infinity. The metric 
outside the shell can be written as 

ds 2 = — ^1 — — ^ dt 2 + ^1 -^jj dr 2 + r 2 dYl 2 (1) 

The interior of the shell is a flat spacetime due to the 
Birkhoff theorem 


ds 2 =-dT 2 + dr 2 + r 2 dfl 2 . (2) 


The time coordinates of the two regions are related with 
the proper time inside the shell as 


dT 

dr 


\J l + R t 2 , 


dt _ \/B + R t ' 2 
dr B 


where B = 1 — R s /R and R T = From here we get 



( 4 ) 


An action of the massless scalar field propagating in the 
background of the collapsing shell can be written as 

S = J d^xy/^g^g^d^dvcf) (5) 
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where <j> is a scalar field, which we can expand in terms 
of the modes as 




In the interior of shell, the action takes the form 

(M) 2 


r rR(t) 

Si n = 27t J dt J drr 2 


Ttidr&Y 


( 6 ) 


(7) 


Similarly, outside of the shell it becomes 


Sout =2ir dt drr 2 

J J R(t) 


■££ + (i-£)«*>’ 


( 8 ) 

The classical equation of motion for this collapsing shell 
near the horizon can be written as [3j 


Rt=±B\ 1- 


BR A 


(9) 


where h is a constant. Using Eq. Q and Eq. Q, we get 
T t = bJi + (1-B)^ (10) 


When the shell is approaching its own Schwarzschild ra- 

n 

comes 


dius, R —> R s , then T t —> 0, hence the total action be- 


r R s 


S 


27T J dty-g J drr 2 (dt<j)) 2 + 
f™ drr 2 (l-^) 


( 11 ) 


which in terms of the modes gives 

S = [ dt f — —— —j—Akk' , k + xQfc-Bfcfc'Qfc' 
/ V 2 B dt dt 2 


with 


( 12 ) 


A kk >=4n drr 2 f k (r)f k >(r) (13) 

Jo 

B kk ' = 4ti- y drr 2 ^1 - f k {r)f k >(r). (14) 

Matrices A kk ' and B kk > are independent of R(t). From 
the action (121, we can find the corresponding Hamil¬ 


tonian and write down the Schrodinger equation Hip = 
idip/dt as 


R \ \ 1 

1 J ) kk 'B k ' + —a k B kk ia k i 


where the momentum is defined as 

d 


n fc = —i 


da k 


dy 

(15) 

(16) 


Since matrices A and B are symmetric and real, the prin¬ 
cipal axis theorem implies that both can be diagonal¬ 
ized simultaneously with respective eigenvalues a and /3. 
One can then write the Schrodinger equation in terms 
of eigenmodes y (which are linear combinations of the 
original modes a) as 


I* 2 

( RJ 2a dy 2 + 2® V 


Defining 


V= I d*(l-^ 


ip(y,t) = i 


Rs 


■ 9ip{y,t) 


dt 


(17) 

(18) 


one can rewrite Eq. in form similar to the harmonic 
oscillator equation as 


1 d 2 

2 a dy 2 




ip(y,r]) = i 


dip{y,y) 

drj 


where 


*<■» = (?) 1 - *S- 


B 


(19) 


( 20 ) 


Comparing Eq. (18) with Eq.(10), it is easy to see that i] 


is basically the time coordinate T for an observer inside 
the shell in the near horizon limit. The exact solution to 
this equation is :B| 


tp(y,v) = e l5iv ' > 


r a 

1 

4 

[. (Or, , 

ay 2 ] 

-7 t9 2 . 

exp 
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where 9 is the solution of the differential equation 

1 


9. 


vv+u 


with initial conditions 

0 ( 0 ) = 


1 

V^o 


,0„(O)=O. 


( 21 ) 


( 22 ) 


(23) 


Since the background spacetime is time dependent, we 
make a distinction between the initial frequency uiq at 
which the mode is created from the vacuum, and the 
final frequency at some later time t defined as 


oj — uj^e^^ 2 , 


(24) 


where we used Eq. (201 and fact that from Eq. ([9]) in the 
near horizon limit we have B = e -t , where the time t 
is expressed in units of Rs- However this frequency u> is 
defined for the time parameter r), so to convert in to the 
the frequency defined for an asymptotic observer, f2(f), 
we have to perform a transform 


o d V- -t- -i/2 

=—uj = e uj = LUoe 1 . 
dt 


(25) 


The wave function ip(y,t) contains information about 
the modes/particles excited in the spacetime in terms 
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of their frequencies at the final moment t. We want to 
construct density matrix of the system so we need to 
expand the wavefunction in terms of a complete basis. 
We will use the simple harmonic oscillator (SHM) basis 
C n(y)- 


V>( V,t ) =^ 2 c n (t)^ n (y) ( 26 ) 

n 

The number of states in this basis is infinite so the size 
of the density matrix will be infinite too. However one 
can see that the probability of exciting higher n states 
decreases rapidly as n increases. Therefore one can eas¬ 
ily identify trends even by considering finite (but large 
enough) n. The coefficients c n (t) can be written as 

Cn(t) = j dyCn*(y)^(y,t). ( 27 ) 


The probability of finding a particle in a particular state 
n is given by | c n (t) | 2 . The coefficients c n can be ex¬ 
plicitly found as (see supplemental material) 


c n (t) 


(—l) TI / 2 e i “ l~2 ( 

(HeV 2 ) 1/4 V P V 


2\" /2 (n — 1)!! 

P J \fn\ 


( 28 ) 


where P is given by 


P = 1 - 


^7 , j_ 

He* V 0 9 2 


( 29 ) 


In order to find c n we need to solve for 6. The simplest 
analytic method is given in [2T1 . 6 and 9 V can be found 
in terms of 77 and £ as 


9 = ^^e+x 2 - 

= (££rj + XXl j)- 

Wo P 


( 30 ) 

(31) 


77 and £ and their derivatives can be written in terms of 
Bessel’s function as 


7 TU , 


i = — [Yo( 2 w 0 )</i(it) - J 0 (2ujo)Y 1 (u)\ ( 32 ) 


TTU . 


x = — [W (2w 0 )Ji(u) - J 1 (2w 0 )yi(u)] (33) 


£»j = -7rwp[y'o(2w 0 )Jo(u) - J 0 {2u} 0 )Y 0 (u)] (34) 


Xv = - 7 rWo[Fi( 2 o; 0 )Jo(u) - J 1 {2u 0 )Y 0 {u)\ ( 35 ) 

where u = 2 ujq^/1 — 77 . 

The occupation number at eigenfrequency H is given 
by the expectation value 

N(t,n) = J2 n K\ 2 . ( 36 ) 


The process of the gravitational collapse takes infinite 
time for an outside observer, however, radiation is pretty 
close to Planckian when the collapsing shell approaches 
its own Schwarzschild radius (see appendix). Since we 
are already working in a near-horizon approximation, if 
we plot N(t, O) for some fixed late t, the spectrum will 
resemble the thermal Hawking distribution [B]. [An im¬ 
portant thing to notice here is that we obtain an ap¬ 
proximately Planckian distribution though we did not 
trace out any degrees of freedom. When calculating radi¬ 
ation from a pre-existing horizon, one gets the Planckian 
distribution of emitted particles only when the infalling 
modes are traced out. In contrast, in the foliation of 
an asymptotic observer all the excited modes remain in 
the causal contact with an observer (even the modes in¬ 
side the shell), and we are not allowed to trace out any 
modes.] While the distribution (36) keeps track only of 


the diagonal terms, we are here interested in correlations 
between the emitted quanta, which is contained not in 
the diagonal spectrum, but actually in the total density 
matrix for the system. 

Density Matrix. Knowing the expansion coefficients 
c n explicitly, we can construct the density matrix. The 
density matrix is defined as 


p=Yl 1 vo w>i 

In our basis it can be re written as 

P — ^ ^ Cmn |Cm) (C n 


(37) 


( 38 ) 


where c mn = c m c n . Original Hawking radiation den¬ 
sity matrix, p^, contains only the diagonal elements c nn , 
while the cross-terms c mn for m ^ n are absent. The 
off-diagonal terms represent interactions and correlations 
between the states. The rationale behind neglecting the 
cross-terms is that these correlations are usually higher 
order effects and will not affect the Hawking’s result in 
the first order. However, as argued recently in [HI: (see 
also [2111), during the process of Hawking radiation, the 
correlations may start off very small, but gradually grow 
as the process continues. It may happen at the end that 
these off-diagonal terms can modify the Hawking density 
matrix significantly enough to yield a pure sate. The 
time-dependent functional Schrodinger formalism is es¬ 
pecially convenient to test this proposal since it gives 
us the time evolution of the system. In Fig. [T] we plot 
some terms (both diagonal and off-diagonal) in the den¬ 
sity matrix. We plot their time evolution with the fixed 
frequency H. We chose a small value of fi to clearly 
demonstrate the effect (which is more pronounced for 
small frequencies). However we should note that Q is an 
exponentially redshifted frequency, and the frequency of 
that mode at the time of creation was much higher. We 
took absolute values of the off-diagonal c mn because they 
can be imaginary. All the units are dimensionless. Di¬ 
mensionless frequency is given as H Rs , while dimension¬ 
less time is given as t/Rs■ From the plot one can see that 
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0 = 0.01, n = 151 t = 5, n = 101 




FIG. 1: Elements of the density matrix c m „ and Tr(p) as 
a function of time at = 0.01, where an index n labeling 
the modes goes up to n = 151. [Note that this is a final 
exponentially redshifted frequency, which was much higher 
at the time of the mode creation.) As time increases, the 
magnitude of coo decreases, Tr(p) remains unity, and all other 
Cmn increase, reach the maximum values then decrease before 
reaching some form of a plateau. This implies that small 
correlations between the modes become as important as the 
diagonal terms. 


the coefficient coo is initially almost exactly one, but then 
it decreases with time. The higher terms start small but 
then they increase with time, reach their maximum value 
and then they decrease before reaching some form of a 
plateau. This is expected because the system starts in the 
ground state. As time progress more modes are excited 
and higher order terms increase in magnitude. This in¬ 
crease of higher order terms can not proceed indefinitely 
if unitarity is preserved, i.e. any increase must be paid by 
a decrease somewhere else. On the same plot, we show 
the trace of the density matrix Tr(p) as a check. The 
trace must remain unity at all times to preserve proba¬ 
bilities. However, we can numerically take into account 
only a finite number of modes. Therefore, at some late 
time, the trace will start decreasing on the graph since 
higher modes which have not been included in numerics 
will become important. The more modes we include, the 
longer the trace will remain unity. In the supplemental 
material we proved that if one takes n —> oo, then Tr(p ) 
always remains unity. Hence we plotted the graph only 
up to the time when Tr{p ) remains one. 

What is more important is that the magnitudes of 
the off-diagonal terms also increase with time. This im¬ 
plies that correlations among the created particles in¬ 
crease with time up to the point when even higher or¬ 
ders terms start increasing. Since there are progressively 
more cross-terms than the diagonal terms, their cumula¬ 
tive contribution to the total density matrix simply can 
not be neglected. In Fig. [2j we plotted c mn and Tr(p) as 
a function of O at a constant time. Tr(p) remains one 
for all frequencies. The lowest term coo increases with 
H, but all other terms decrease. This means that the 


FIG. 2: Cross-terms c„m and Tr(p) as a function of f 1 at 
fixed time t = 5. As Q. increases all c m „ decrease, but coo 
increases. 


lowest diagonal term dominates and correlations are not 
that important at high frequencies. Information content 


Q = 0.0005, n = 401 



FIG. 3: ph is the diagonal Ha wki ng density matrix, p is the to¬ 
tal density matrix as in Eq. (381. We plot Tr(p 2 ) and Tr(p^) 
as functions of time at a fixed final frequency = 0.0005. 
The magnitude of 7Y(p|) decreases with time which means 
that the system is losing information by going from a pure 
to a mixed state. However Tr(p 2 ) remains unity at all times, 
which means that the state remains pure. This implies that 
the information of the system is conserved if cross-correlations 
are accounted for. 


in the system is usually given in terms of a trace of the 
squared density matrix. If the trace of the squared den¬ 
sity matrix is one, then the state is pure, while the zero 
trace corresponds to a maximally mixed state. In Fig. [3j 
we plot the traces of squares of two density matrices as 
functions of time for a fixed frequency. One is the Hawk¬ 
ing density matrix ph which contains only the diagonal 
terms c nn and neglects correlations. The other one is the 
total density matrix p defined in Eq. (381 which contains 
all the elements, including the off-diagonal correlations. 
As expected, Tr{pj l ) decreases as time progress which 
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means that the system is going from a pure state to a 
mixed thermal state. This is often labeled as the infor¬ 
mation loss in the process of Hawking radiation. How¬ 
ever, if the plot the total Tr(p 2 ) we see that it always 
remains unity, which means that the state always remain 
pure during the evolution and information does not get 
lost. This clearly tells us that correlations between the 
excited modes are very important, and if one takes them 
into account the information in the system remains in¬ 
tact. In Fig. [4j we plot Tr(p 2 ) and Tr(p\) as a function 


t = 5, n = 51 



FIG. 4: Tr(p 2 ) and Tr(pf l ) as a function of fi at i = 5. Again 
Tr(p 2 ) remains one at all frequencies but Tr(p\) differs from 
unity at low frequencies which means it does not account for 
the full information of the system at low frequencies. 


of fl at a fixed time. As expected, Tr(p 2 ) remains one at 
all frequencies, but TV(/5^) differs from unity at low fre¬ 
quencies. This implies that pu gives a good description 
of the system at high frequencies, but it fails to do so at 
low frequencies. 

Conclusions. In conclusions, we showed by explicit 
calculations that radiation coming from a collapsing ob¬ 
ject is manifestly unitary in a space-time foliation that 
is motivated by an asymptotic observer. Hawking’s ther¬ 
mal density matrix is diagonal and inevitably leads to in¬ 
formation loss. However, when we take the off-diagonal 
correlation terms into account, the density matrix de¬ 
scribes a pure state at all times. This result agrees well 
with ] 22 j , where it was shown at that at late enough time 
all the information in the system is contained in correla¬ 
tions between the small subsystems (in this case emitted 
particles). For unitarity to be manifest to an observer, he 
has to be able to observe the total density matrix, i.e. all 
the created modes and correlations between them. In the 
foliation of an asymptotic observer all the excited inodes 
remain in the causal contact with an observer (even the 
modes inside the collapsing shell since the collapsing ob¬ 
ject must disappear in finite time as seen by this ob¬ 
server), and we are not allowed to trace out any modes. 
This is in contrast with the usual procedure when Hawk¬ 
ing radiation is calculated from a pre-existing horizon or 
in the limit when the horizon is formed at infinite future 


where the ingoing modes become causally separated from 
the outside observer. 

If some of the modes are lost into the singularity (e.g. 
in the standard calculations like in E3)> then the incom¬ 
plete density matrix may not look like that of the pure 
state. However, there are strong indications that quan¬ 
tum effects should be able to rid gravity of singularities, 
just as it was the case with the singular Coulomb poten¬ 
tial (see for example 123]). 

Finally, our analysis was done for a static outside ob¬ 
server, however, it will be very important to learn what 
an infalling observer would see during the collapse in or¬ 
der to settle down the question of information loss. 
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Appendix A: Number of particles produced as a 
function of time 

We work in a space-time foliation of an asymptotic ob¬ 
server. This observer measures frequency Q according 
to his time parameter, t. From the original Schrodinger 
equation (Eq.17) expressed in the asymptotic time co¬ 
ordinate, f, we see that the mass term of the harmonic 
oscillator is a(t) which is related to a constant a as 

a(t) = . ^ = «e‘. (Al) 

V 1 R ) 

To go to Eq.(19), we divide (Eq.17) with the factor B = 
1 — R s /R. From Eq.(9) we have B = e - *, where the 
time t is expressed in units of Rs■ This will rescale our 
time parameter from t to 77 . We then proceed by solving 
Eq.(19), and the solution is given by Eq.(21). A solution 
to Eq.(19) must also be a solution to Eq.(17) with an 
appropriate rescaling. 

To study the particle content in the system, we need 
to expand the wavefunction in Eq.(19) in terms of a com¬ 
plete basis. We use the simple harmonic oscillator (SHM) 
basis C n(y)- 


i>{y,t) =^2c n (t)Cn(y) (A 2 ) 

n 

The coefficients c n (t) can be written as the overlap be¬ 
tween the wavefunction and the basis states 

c n (t) = J dyC n *(y)ip(y,t). (A3) 

The simple harmonic oscillator basis functions in terms 
of the frequency f l measured by an asymptotic observer 
are 


/a(t)fl \ 1 ^ 4 e- a ^) n v 2 / 2 , _ 

Uv) = ( - 11 - H n (^amy) (A4) 

\ it J V2 "n! 
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where H n are the Hermite polynomials. A mode is cre¬ 
ated from vacuum at its initial frequency ujq = (as 
in Eq. (20)). So in vacuum, Q = uiq. Since the back¬ 
ground is time-dependent, this frequency will evolve in 
time according to Eq.(25). Then, Eq.(21) together with 
Eq.(26) gives 



To find /„ consider the corresponding integral over the 
generating function for the Hermite polynomials 


,J(z) = J dxe- Px2/2 e~ z2+2zx 


0 -z 2 (l-2/P) 


(A7) 


Since 


e -^ +2 z, = y- z Hn{x) 

' n\ 

n—0 


(AS) 


dxe Px2/2 H n (x) = — J(z) 


z=o 


(A9) 


Therefore 

I o_ / 9 \ 

In = \lp\ 1 -p) Hn{0) ■ (A10) 

Since 

H n { 0) = (All) 

Vn! 

and H n ( 0) = 0 for odd n, we find the coefficient c n for 
even values of n, 


_(-l) n / 2 e iS l~2 / 2 \ " /2 (n — 1)!! 

Cn “ (ne‘0 2 ) 1 /* Vp(}-p) Vn! 

For odd n, c n = 0. 


(A12) 


Appendix B: Planckian spectrum 

In his original calculations, Hawking used the Bo- 
goliubov transformation between the initial (Minkowski) 
vacuum and final (Schwarzschild) vacuum at the end of 
the collapse. The vacuum mismatch gives the thermal 


spectrum of particles. In this picture, there is a negative 
energy flux toward the center of a black hole and positive 
energy flux toward infinity, and thus a black hole loses its 
mass. Note that the existence of the horizon is necessary 
for this, since it is not possible to have a macroscopic 
negative energy flux without the horizon (the fact that 
the time-like Killing vector becomes space-like within the 
horizon is responsible for this). Since an outside observer, 
which is the most relevant observer for the question of the 
information loss, never sees the formation of the horizon, 
it is not clear how this picture works for him. 

Here we plot the distribution of produced particles 
N(t, H) = J2 n n \ c n\ 2 f° r several different times, and com¬ 
pare it with the Planckian distribution -/Vpi anc k(fi) = 
l/(e n / T — 1), where T is the temperature. 


N 



FIG. 5: Particle occupation number N(Q) for several differ¬ 
ent times. As time increases, the distribution approaches the 
Planckian distribution Api an ck(H) (in particular, N(Q) be¬ 
comes divergent at SI —» 0). 


This is where the advantage of the time dependent 
functional Schrodinger formalism becomes obvious. In 
this picture, it is the time dependent metric during the 
collapse which produces particles. No horizon is needed 
for this. As the collapsing shell approaches its own 
Schwarzschild radius, this radiation becomes more and 
more Planckian. Finally, at t -A oo, when the horizon 
is formed it becomes completely Planckian and we re¬ 
cover Hawkings result (his calculations were performed 
in t —> oo). 

Note that a static outside observer will never witness 
formation of the horizon since the collapsing object has 
only finite mass. An outside observer will observe the 
collapsing object slowly getting converted into Hawking 
radiation before horizon is formed. For him, no horizon 
nor singularity ever forms. Still, if the radiation com¬ 
ing out from the collapsing object is thermal, and since 
there is a region inaccessible to an outside observer in 
Schwarzschild coordinates, this would lead to informa¬ 
tion loss. However, if radiation is not thermal, then in¬ 
formation from the part of the space inaccessible to an 
outside observer could leak out. 
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We now emphasis that the Planckian spectrum of pro¬ 
duced particles is not equivalent to a thermal spectrum. 
For a strictly thermal spectrum there should be no corre¬ 
lations between the produced particles. The correspond¬ 
ing density matrix should only have non-zero entries on 
its diagonal. In contrast, if subtle correlations exist, then 
the distribution might be Planckian, but the density ma¬ 
trix will have non-diagonal entries. If the state is pure, 
then the trace of the squared density matrix will be unity. 

It is crucial for an observer to be able to see the whole 
density matrix, i.e. that he is able to measure all the 
modes and their correlations at least in principle. In the 
foliation of an asymptotic observer that we worked in, all 
the excited modes remain in the causal contact with an 
observer (even the modes inside the collapsing shell), and 
we are not allowed to disregard any modes. In contrast, 
when calculating radiation from a pre-existing horizon, 
one has to trace out the infalling modes which eventu¬ 
ally get lost into the singularity. An important thing to 
notice here is that we obtained an approximately Planck¬ 
ian distribution though we did not trace out any degrees 
of freedom, which is again in contrast with calculations 
in the presence of the pre-existing horizon where ther¬ 
mal spectrum is obtained after tracing out the infalling 
modes. Even in that case, thermality might be only ap¬ 
parent for some limited time. Singularity at the center 
is a classical result. Most likely, singularity at the center 
can be cured by quantization [23], just like we cured the 
atom of the classical 1 jr singularity of the electrostatic 
potential. If there is no singularity at the center, and 
instead of it we find only a region of very strong but fi¬ 
nite gravitational fields, then no modes are lost for an 
outside observer forever, and he can again measure the 
whole density matrix. 


Appendix C: Trace of the density matrix 


Here we will prove that the trace of the density matrix 
will add up to unity if n —> oo. Let’s take 


5 



(Cl) 


Now, 


Tr(p)= Y, 


n=even 

2 




VQePPlPl ^ n\\ 

v I I n=even 

2 1 

VsW\p\ 

2 1 




L i 1 p i 


(C2) 


Inserting the expressions for P and little algebra gives 
Tr{p) = 1- 

It is instructive to check also how the number of excited 
modes n that we include in the calculations influences the 
behavior of Tr(p), Tr{p\) and Trip 2 ) with time. The 
trace of the original density matrix Tr(p) must remain 
unity at all times to preserve probabilities. However, we 
can numerically take into account only a finite number of 
modes. Therefore, at some late time, the trace will start 
decreasing on the graph since higher modes which have 
not been included in numerics will become important. 
The more modes we include, the longer the trace will re¬ 
main unity. We can also see in Fig. [6] that Trip 2 ) remains 
unity as long as Trip) is unity. However, the Hawking 
Trip\) drops from unity much earlier, and several curves 
with different n lie on top of each other. Since ph has 
only diagonal terms, this tells us that the off-diagonal 
correlation terms indeed play crucial role in preserving 
purity of the state. 



FIG. 6: Tripi) an d Trip 2 ) are plotted for different values 
of n. Trip 2 ) remains unity as long as Trip) is unity. How¬ 
ever, the Hawking Ti'ip \) drops from unity much earlier, and 
increasing n does not affect it much (several curves with dif¬ 
ferent n lie on top of each other). 
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